home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / cpptrs.z / cpptrs
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCPPPPPPPPTTTTRRRRSSSS((((3333FFFF))))                                                          CCCCPPPPPPPPTTTTRRRRSSSS((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CPPTRS - solve a system of linear equations A*X = B with a Hermitian
  10.      positive definite matrix A in packed storage using the Cholesky
  11.      factorization A = U**H*U or A = L*L**H computed by CPPTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE CPPTRS( UPLO, N, NRHS, AP, B, LDB, INFO )
  15.  
  16.          CHARACTER      UPLO
  17.  
  18.          INTEGER        INFO, LDB, N, NRHS
  19.  
  20.          COMPLEX        AP( * ), B( LDB, * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      CPPTRS solves a system of linear equations A*X = B with a Hermitian
  24.      positive definite matrix A in packed storage using the Cholesky
  25.      factorization A = U**H*U or A = L*L**H computed by CPPTRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              = 'U':  Upper triangle of A is stored;
  31.              = 'L':  Lower triangle of A is stored.
  32.  
  33.      N       (input) INTEGER
  34.              The order of the matrix A.  N >= 0.
  35.  
  36.      NRHS    (input) INTEGER
  37.              The number of right hand sides, i.e., the number of columns of
  38.              the matrix B.  NRHS >= 0.
  39.  
  40.      AP      (input) COMPLEX array, dimension (N*(N+1)/2)
  41.              The triangular factor U or L from the Cholesky factorization A =
  42.              U**H*U or A = L*L**H, packed columnwise in a linear array.  The
  43.              j-th column of U or L is stored in the array AP as follows:  if
  44.              UPLO = 'U', AP(i + (j-1)*j/2) = U(i,j) for 1<=i<=j; if UPLO =
  45.              'L', AP(i + (j-1)*(2n-j)/2) = L(i,j) for j<=i<=n.
  46.  
  47.      B       (input/output) COMPLEX array, dimension (LDB,NRHS)
  48.              On entry, the right hand side matrix B.  On exit, the solution
  49.              matrix X.
  50.  
  51.      LDB     (input) INTEGER
  52.              The leading dimension of the array B.  LDB >= max(1,N).
  53.  
  54.      INFO    (output) INTEGER
  55.              = 0:  successful exit
  56.              < 0:  if INFO = -i, the i-th argument had an illegal value
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.